Boundary contributions to the hypervirial theorem 
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It is shown that under certain boundary conditions the virial theorem has to be modified. We 
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harmonic oscillator with point interaction have been studied in the light of this generalization of the 
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I. INTRODUCTION 

The virial theorem in classical and quantum mechanics 
has been known for a very long time. It appears for the 
first time in the classical works of Clausius in 1870 [l| 
in the context of statistical mechanics. See also 0] for a 
historical account. Already in 1930, Fock Q derived the 
quantum mechanical version of the theorem that relates 
the expectation value of the kinetic energy in an eigen- 
state V'n of the Hamiltonian to that of the Clausius virial 
function. More precisely 



(1) 



The theorem was initially derived for eigenstates of non 
relativistic Hamiltonian in flat space with potentials de- 
pending only on the position, but it can be generalized in 
many different ways: for relativistic systems, for a par- 
ticle moving in an electromagnetic field, to many body 
systems, mixed states and for compact space (H-pjJ to 
mention a few. Also recently it has been formulated a lo- 
cal version of the virial theorem for systems of fermions, 
see refs. [12] and [l3j . 

Another kind of generalization is the hypervirial theo- 
rem [l~i| that works, in principle, for any operator G and 
can be written 

(V> n |[ff,G]^} = 0. 
The virial theorem in dimension d is recovered when 

i d 

In this paper we want to discuss the domain of validity 
of these theorems and, for the cases in which they fail, 
we show how to modify them by adding the appropriate 



term related to the boundary conditions of the quantum 
system. In this way we extend the applicability of the 
virial and hypervirial theorem to these systems. 

In concrete terms, assume that the Hamiltonian is de- 
fined over a domain D(H). ip n necessarily belongs to 
D(H), but it may happen that Gip n £ D(H). In this 
case the hypervirial theorem has to be modified by the 
addition of an extra term A: 



with 



{^ n \[H,G]^ n ) = A, 



A = (ip n \(H* - E n )Gip n ), 



(2) 



where E n is the energy eigenvalue of ip n and is the 
closed extension of H to a domain that includes Gip n . Of 
course, the new term cancels if = H i.e. when Gip n € 
D(H), but if that is not the case it is different from zero 
in general. We will see that A can be computed as the 
integration of a pure derivative term and, therefore, it 
depends only on the value of the wave function and its 
derivatives at the boundary. This fact is reminiscent of 
the analogous objects in quantum field theory like the 
chiral anomaly [15]. We would like to remark that an 
additional term for the virial (or hypervirial) theorem in 
classical mechanics may appear also when the orbits are 
not bounded in phase space. 

This paper is organized in the following fashion: in 
the next section we discuss a simple example to motivate 
our study and obtain the generalization of the hyper- 
virial theorem. Then we give two examples of systems 
where this generalization is needed, in section III we dis- 
cuss the Hydrogen atom problem with point interaction 
in three dimensions while the section IV is devoted to 
the isotropic harmonic oscillator. Finally we sketch our 
conclusions in section V. 
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II. GENERALIZED HYPERVIRIAL THEOREM 

Before discussing the required generalization of the hy- 
pervirial theorem we would like to motivate our study by 



2 



considering a very simple example. 

Consider a free particle in one dimension restricted to 
move in [0, oo) and subject to Robin boundary condi- 
tions, i.e. ^'(0) + aip(0) = with a > 0. In this case the 
free Hamiltonian has a single eigenfunction 



with eigenvalue 1 



En 



2a e~ 



2m 



If we use (Q} to compute the expectation value of the 
kinetic energy in this state we obtain that it vanishes, 
which is in contradiction with the real result 



(MTih) = e . 



(3) 



The reason for this apparent contradiction is the fact 
that the domain of the Hamiltonian is not preserved by 
the generator of the scale transformation 



r 1 1 
G= —h Xp -2 



(4) 



and the virial theorem has to be modified. 

To understand the origin of this modification let us 
consider in detail the derivation of the virial theorem. 
We assume a quantum system with infinite dimensional 
state space and a self-adjoint Hamiltonian with domain 
D(H). If Gtp n G D(H) we have 



(*l> n \[H, G]lp r , 



(^ n \HGlP r , 

(H<p n \G^ 



(lPn\GHA, 
{^ n \GH^ r , 



0(5) 



where in the second line we have used that H is self- 
adjoint and that ip n is its eigenstate. 

The problems appear if Gip n g" D(H). In this case the 
expressions in the first line of (J5J) do not make sense unless 
we extend H to a larger domain that includes Gip n . In 
principle it can be done in an arbitrary way but in this 
case there is a well defined prescription to obtain the 
required extension. First we restrict H to an appropriate 
dense subspace S, typically to functions whose support 
does not include the boundary. If we define H — H\s, the 
extension H„ is the adjoint of the previous restriction, i.e. 

= H + . We will assume that the restriction is closed 
(which could always be achieved by taking the closure of 
its graph) and then we know = H. 



1 At first sight it may seem strange that T = p 1 /(2m) has a neg- 
ative eigenvalue. Note, however, that p is not symmetric with 
Robin boundary conditions and in fact one has 

(m) = £-«tfW> +^(o)V(o» = ^-ww) - a\m\ 2 ), 

2m 2m 

and the second term on the right provides a negative contribu- 
tion. 



With the previous machinery we can rewrite ([5]) in a 
meaningful way 



(1> n \[H.,G\ih 



(>>p n \H*GTp n ) - {^ n \GH^ n 

(lpn\(H* — E n )Glp n ) = A, 



and obtain the generalized virial theorem as announced 
in the introduction. Looking at the last line of this ex- 
pression one could be tempted to take the adjoint of the 
operator H* and make it act in the left side of the scalar 
product. However, this is not possible in general as the 
adjoint of is H with a domain smaller than H that 
may not contain ij) n . Of course, if it so happens that 
ip n £ S the additional term cancels. 

We would like to show now that under certain assump- 
tions (that we actually meet in our examples) the extra 
term is a boundary term, i.e. it depends only on the value 
of the functions and their derivatives at the boundary. 

First of all, we rewrite the extra term in a way that 
is appropriate for wave functions that are not necessarily 
eigenvectors of H, 

A = (i/> n \H.Ghl> n ) - (H^ n \G^ n ). 

Then we will show that the value of A does not change 
if we replace if) n by ip n + x with x € S and G\ € S. The 
difference reads 



A A = (il> n \H*G X ) 
+ (xl#*Gx)~ 



- (Hi; n \Gx)- 

- (H x \G^ n )- 
(H X \G X ). 



(6) 



Now, given that G\ € E>(H), we can replace iJ» by H 
in the first term of this expression and, if we take its 
adjoint, it cancels the second term. For the very same 
reason the last two terms cancel. Also in the third term 
we can take the adjoint of if*, because x lS m t ne domain 
of = H, and therefore it cancels the fourth term and 
A A vanishes. 

Then if, as we supposed before, the elements of S 
are wave functions whose support does not contain the 
boundary we deduce that the extra term A is invariant 
under deformations of that do not affect the values 
of the function and its derivatives at the boundary or, in 
other words, the extra term depends only on the latter. 

As an illustration to these results and before discussing 
the examples of physical interest, we will compute the 
new term for the simple system introduced at the begin- 
ning of this section. One immediately gets 



A 



2m 



2m 



(Chh(x)^(x) - 1>S(x)%Ghl> (x)) dx 



d x (Gip (x)d x %{x)- 
-^Q(x)d x Gip {x))dx 



(7) 



which is a pure boundary term (for any operator G) as 
stated. 

In our particular case it amounts to 

A = £ (^o(O)^(O)* - Vo(0)Vo(0)*) = 



3 



fi. 2 a^ 2 h 2 a 2 
2m m 



(8) 



Now taking into account that [G, H] — 2T and plugging 
these expressions into the generalized hypervirial theo- 
rem J2|) we get 

2(V |7V>o) = 



from which © follows. 

We have seen that if we modify appropriately the hy- 
pervirial theorem it can be applied to situations in which 
the action of the operator G on the energy eigenstate puts 
it out of the domain of the Hamiltonian. In the following 
sections we will apply this idea to the Coulomb problem 
and later on to the isotropic harmonic oscillator. 



III. COULOMB POTENTIAL. 

In this section we shall discuss the virial theorem for 
the Coulomb problem in three dimensions. We shall see 
that with the chosen boundary conditions at the origin 
the extra term is required to cancel the divergence that 
appears in the application of the virial theorem. 

Consider the Hamiltonian 



H 



P 

2m 



(9) 



If we separate the angular variables from the radial one 
we obtain the Hamiltonian for the latter 



Hi = 



2m 



1<P_ 
r dr 2 



1(1 + 1) 



(10) 



Here 



that depends on the angular momentum number I. 
£ = 2mk/h 2 . 

We can simplify the Hamiltonian by performing a sim- 
ilarity transformation <f> — rip. Therefore the new wave 
functions <fi are square integrable in [0, oo) with respect 
to the Lebesgue measure. The new Hamiltonian is then 



H, 



2m 



dr 2 r 2 r 



(11) 



and it is symmetric with respect to the Lebesgue mea- 
sure in [0,oo) when acting on normalizable wave func- 
tions whose support does not contain the origin. 

In order to have a self-adjoint Hamiltonian we have 
to look for extensions of Hi with such a property. This 
issue has been analyzed in detail in ([16]). The result is 
that if I > 1 there is a single self-adjoint extension of Hi 
whose domain are the wave functions that vanish at 0, 
while Hq has an infinite number of self-adjoint extensions 
characterized by a single parameter a. If we call H a such 
a self-adjoint extension, its domain is given by 



D{H a ) = {(/) <= D(H+)\ as r -> 0, 

<Kr) = 0(O)(l-£rln(|£|r) 



where the little-o notation has been used. The case of 
wave functions that vanish at the origin (the standard 
boundary conditions for the hydrogen atom) is recovered 
when a — > — oo and (j>(0) — > 0, while the product (f>(0)a 
remains finite. 

The point-spectrum of H a is obtained from the solu- 
tions in A to the equation 



1 



F c (\) = *(1 - A) - In |A| + — + 2 7 - 1 = 



(12) 



such that £/A > 0. Here 7 is the Euler's constant and $ 
is the digamma function, i.e. the logarithmic derivative 
of the gamma function. For a solution A„ of the previous 
equation the eigenvalue of H a is 



1 



2h 2 XI 



In fig. 1 we represent graphically the solutions of (TT2"j) for 
A, as the intersection of the branches of the curve with 
the line of constant value a/£. In the negative side F(X) 
has an asymptotic value 27 — 1 s» 0.1544. 




FIG. 1: In the figure we represent A in the horizontal axis 
and we plot Fc(X). The point-spectrum of H a is obtained 
from the values of A at the intersection of the graph of Fc 
with the horizontal line at a/f. Only those values of A with 
the same sign that £ are allowed. 

Note that even in the presence of a repulsive Coulomb 
potential, k < 0, if we take a/£ < 27 — 1 there is a single 
bounded state associated to the negative solution for A 
of the equation Fc(X) = a/£. 

The eigenfunction associated to the eigenvalue E n is 
given by the Whittaker function [l7j conveniently nor- 
malized 

<t>n{r) =N n W XnA/2 (£r/X n ). 
With the normalization constant given by 

£|r(-A„)| 2 



o(r))}. 



N z = - 
" 2A„*'(-A„) + 2 - Xn 1 ' 

where the prime denotes derivative. 



(13) 
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Now we want to apply the generalized form of the 
virial theorem to this system. We take the operator 
G = — 2^xp — |, then the modified hypervirial theorem 
reduces to 

2(i/) n \Tip n ) + (ip n \Vip n ) = A. 

Using the energy eigenvalue and integrating out the angu- 
lar coordinates we can write for the radial wave function 



A 



r 



2E n 



(14) 



In this expression we can immediately see the necessity of 
an additional contribution to the virial theorem: the ex- 
pectation value of the potential term diverges logarithmi- 
cally at the origin (recall that we must use the Lebesgue 
measure in [0, oo) ) while the right hand side is finite. In 
fact, we shall see that an analogous divergence appears 
in the extra contribution so that both of them cancel to 
yield the correct finite result. 

We will sketch how to verify it. We first introduce 
a cut-off e that removes the origin from the integration 
region. The extra term only depends on the boundary 
conditions of the wave function. We can take it from ([7]) 
and in its regularized version reads 

A =-!H^(0)| 2 (£ln(|e|e)+£- « + •••), 
2m 

where the dots stand for terms that vanish when e — > 0. 
We now compute the potential term 

k f°° 1 

(ct>n\-<l>n)e = N 2 k / -W Xn >1/2 (£r / A„) 2 dr. 

To evaluate the expectation value we use 



-W XA/2 {xf dx = 

' (Ine - A#'(l - A) + f(l - A) + 2 7 + • • • ) . 



r(-A) 

where, as before, the dots represent terms that vanish 
when e — > 0. The integral can be obtained using the 
identity [H| 



-W^(x)W p ^)dx 



1 



(i-p 



[W^{x)W' p! „{x) - W'^(x)W p , a (x)}(15) 



And the series expansion of W M „• [171 ]. 
Putting everything together and given 

W A ,i/a(0) = l/r(l-A) 
and (|12[) we obtain, 



A 



k 



r 



mk 2 1 

"~^aI 



in agreement with (|14[) . 



IV. HARMONIC OSCILLATOR 

In the previous section we saw that the verification of 
the virial theorem were plagued by infinities. This fact 
is not inherent to the generalized virial theorem itself 
but it is rather related to the singularity in the Coulomb 
potential. To illustrate this fact we shall discuss now a 
new instance of the generalized virial theorem in three 
dimensions that is free of these divergences. 

We consider a 3 dimensional harmonic oscillator given 
by the Hamiltonian 



H = 



2m 



1 2 2 

— mu> r 
2 



(16) 



As before, we factorize the angular part and study the 
radial Hamiltonian for the radial part ip(r) of the wave 
function of angular momentum number I: 



H, = 



2m 



r dr 2 



+ , 1 2 2 

+2 m0Jr ■ 



Or acting on <p(r) = rip(r) 
Hi 



2m 



dr 2 



1(1 + 1) 



2 2 

■mw r , 



that is a symmetric operator with respect to the Lebesgue 
measure in [0, oo) when acting on functions whose sup- 
port does not include the point 0. 

In order to obtain a self-adjoint operator we have to 
extend Hi to a larger domain. As it was the case in the 
previous section if I > 1 the only self-adjoint extension 
is to functions that vanish at the origin. In this case 
the extra term cancels and the standard virial theorem 
holds. However, if I — 0, we have a whole family of self- 
adjoint extensions parametrized by /3. If we call H° such 
an extension its domain corresponds to Robin boundary 
conditions, i.e. 



{0e^(i? o +)|0'(O) + 2/^(0) = 0}. (17) 



The normalized eigenfunctions of H@ can be written 
in terms of the Whittaker function as 



<t>E{r) 



N 



W KA/4 (a 2 r 2 ) , 



(18) 



where a = (mu> /K) 1 / 2 , n — and the normalization 
constant N is obtained from the condition 



N 2 [-oo -, 

TT / -W Kil/i (x) 2 dx = 1. 
2a Jo x 



The integral can be performed using (|15|) and finally we 
get 



N 



, _ 2a r(3/4- k)T(1/A-k) 
~ ~*(3/4- k) - #(1/4- «■)' 



5 



Only for specific values of k the eigenstates (TT8l) belong 
to domain (fTT|). demanding ^ g D(H 13 ) we obtain the 
equation for k 



F ff (/s) 



_ r(3/4- k) _ p 



r(l/4-/t) cr ' 

The solutions are represented in fig. 2. 



(19) 



-4 -2 



\ 



FIG. 2: In the figure we represent Fh(k) with k in the 
horizontal axis. The point spectrum of are obtained from 
the values of k at the intersection of the graph of Fh with the 
horizontal line at P/a. The possible values of k are given as 
the intersection of the curve with the horizontal line at j3/a. 

In general the eigenvalues E n = 2n n hu) must be com- 
puted numerically. One particular case in which we have 
analytic solutions is ft = 0, that correspond to Neumann 
boundary conditions. In this case we get E n — nfku + 1/2. 
We also have exact results for j3 = — oo, Dirichlet bound- 
ary conditions and E n = nfrio + 3/2. The latter cor- 
responds to the harmonic oscillator without point in- 
teraction. Note also that for /3/a > r(3/4)/r(l/4) » 
0.337989..., there is an eigenstate with negative energy. 
Of course, as discussed in the footnote of section II, there 
is not a contradiction between the negative eigenvalue 
and the fact that H@ is the sum of squares of opera- 
tors. The point is that with Robin boundary conditions 
some of those operators are not self-adjoint and there is 
not any reason to argue that H@ is semi-positive defi- 
nite. The same applies to the Coulomb potential with 
repulsive interaction that we discussed in the previous 
section. 

The generalized virial theorem, in terms of the radial 
wave function </>„, reads 



{<t>n\T4> n ) - {4>n\V4>r, 



A 



were we have taken G as before. Using the energy eigen- 
value we can rewrite the identity above in the following 
way 



and this is the expression that we will check. 
We first compute the additional term 

a = £-MWn(or = -Si^(°)i 2 

4m 2m 



-N 



1 



2m r(3/4 - k„)T(1/4 - « n ) 



(21) 



In the last equality condition (flUf was used. 

The expectation value of the potential can be written 

ah 2 f°° 

(<P n \Vcfr n ) =N 2 — J W KnA/i (z) 2 dz, (22) 

where the change of variable z = a 2 x 2 was performed. 

Now we can compute the previous integral with the 
help of the identity 

zW K>1/i (z) = W K+ltl/i {z) + 2 K W Ktl/i (z)+ 

+ (3/4 - k)(1/4 - K)W re _ lil/4 (z),(23) 

that we apply to one of the factors of the square inside 
the integral. We perform the resulting integrals using 
(fT5j) and finally we obtain 



{(j>n\V4>n 



<7 2 h 2 K r , 



+N 2 



,irah 2 



4m T(3/4 - K„)r(l/4 - K n ) 



(24) 



2{(j> n \V(l> r , 



A = E n 



(20) 



If we plug this expression and (|2"Tj) into (12^1) we can see 
that the identity holds. 



V. CONCLUSION 

We proved that, when the dilations symmetry does not 
keep invariant the domain of definition of the Hamilto- 
nian, the exact virial theorem has an extra term that 
accounts for this fact. We have shown that this term is 
the integral of a pure derivative and, therefore, its value 
depends only of the behaviour of the wave function at 
the boundary. We have verified this generalization of the 
virial theorem in the cases of Coulomb and harmonic os- 
cillator with point interaction in three dimensions. 

As an extension of our work, it could be interesting to 
study if a similar phenomenon occurs for the local virial 
theorem for fermions with harmonic potential [Hj][l3j], 
and the boundary conditions used in this paper. 
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